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Abstract 

The path-integral calculation for the free energy of a fermion gas in flat spacetime is performed in 
spherical coordinates, and its equivalence with Cartesian coordinates is established. An appropriate 
generalization of spherical harmonics for fermion fields is used. The latter technique has been 
developed by the authors of this paper to perform the calculation of the free energy of a fermion 
gas in the presence of a Schwarzschild background, using a near-horizon expansion approach. The 
robustness of this equivalence in the flat-spacetime limit lends support to the black hole background 
computation leading to black hole thermodynamics. 
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I. INTRODUCTION 



The authors of this paper have recently advanced a near-horizon (NH) expansion tech- 



nique for 
approach 



;he calculation of black-hole entropy in the framework of 't Hooft's brick-wall 
ll. This NH technique has provided a clear origin for the leading divergence of 



the entropy calculation, as well as revealed a NH conformal symmetry that may have an 
important role in the understanding of the Bekenstein-Hawking area law j2-4|. In Ref. ^ an 
expansion in spherical harmonics was used that led to a radial one- dimensional semiclassical 
analysis for the expression of the density of modes. ^ As in 't Hooft's work, this was done 
for spin-zero fields. 

We have recently performed a similar calculation for Dirac fermion fields |^, in both the 
canonical and path-integral frameworks. The NH method in spherical coordinates proved 
to be a critical technique to tackle the technical challenges posed by the fermionic path 
integral in a spherically symmetric black hole background. The Bekenstein-Hawking area 
law was verified within the framework of the brick-wall approach appealing to the relevant 
fermionic spherical harmonics expansion. In this paper we develop for fiat spacetime the 
corresponding expansion in fermionic spherical harmonics, and present further evidence for 
the technical correctness of the subtle fermionic NH approach. While the path-integral 
calculation of the free energy for a free-fermion gas in fiat spacetime is well known using 
Cartesian coordinates jsl, [J, its counterpart in spherical coordinates is not. Through a 
series of steps leading to the free-energy path integral in spherical coordinates, we show 
the equivalence with Cartesian coordinates, and provide a strong check on the correctness 
of the corresponding fermionic path-integral in a black-hole background. In essence, this 
calculation highlights some conceptual and technical properties that are crucial for similar 
computations in curved spacetime, and lends itself to generalizations to higher spin fields. 

Consider the Euclideanized action for free fermion fields in flat spacetime 



Se = J dT J d'x^{Y^d^ + m)ij (1.1) 
(the relevant definitions and conventions are summarized in Appendix |A]) . The partition 



^ The authors of reference |4| have also developed NH methods without the need for an expansion in 
spherical harmonics, and obtained the same results 
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function for this system is 



Z = y D'ip (t, x) (t, x) exp {—Se} ■ 



;i-2) 



:i.3) 



As is well known, the fermionic fields satisfy anti-periodic boundary conditions in imaginary 
time r, 

ip (0, x) = —ip {(3, x) 
^PiO,x) = -^ 

This implies that the frequencies associated with the Euclidean time r are discrete, and of 
the form 

w„ = y + , n G Z. (1.4) 

In flat spacetime, one can always introduce the momentum space Fourier representation for 
the spatial part of the fermion field 

exp (— zcj„r) f (Pp 



The Euclideanized action then becomes 



exp z P -x) ipn ip) ■ 



(1.5) 



:i.6) 



where the four vector P = {un, p)- Therefore, using the well-known result for ant i- commuting 
path integrals, 

j \^da*da^ exp {-a*Mijaj) = det (M) , (1.7) 
the partition function (11. 2p gives 

Z = Cj]det(^7^PM + ^) 

1/2 



n,p 



c 



J]^ det (^Te-Pm + ^) *i6t (-i7^P^ + m) 



n.p 



1.8) 



with C a constant. Here we used 



|detL>| = Vdet D det 



Thus, up to a phase, we get 



Z = \{det y^[{ujl+p^ + m')h] 



:i.9) 



(1.10) 
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where I4 is the 4x4 unit matrix arising from the Dirac algebra. Notice that the determinant 
involves products with respect to the Dirac matrices in addition to its functional nature; 
here, due to the simple factorization above, the matrix part leads to an overall exponent of 
4 for the functional determinant, 

Z = njdetV^K + p^ + m^)]^ 

"'^^ (1.11) 
= JJ det 2 {col + P^ + m^) ■ 

n,p 

From Eq. (11. lip , it is straightforward to get the free energy using F = —\nZ/p. 



II. PATH INTEGRAL IN SPHERICAL COORDINATES 



In spherical coordinates a simple change of coordinates from cartesian to spherical coor- 
dinates changes the Dirac operator in the action to 



= 1%do+ j^rdr+ 



1 



r sin6' 



d^ + m. 



where 



fit } ^fii 



0, Ci ) , i = r,e,, 



(2.1) 



(2.2) 



The operator B, defined above in the original field representation ip (r, x), is not convenient 
for our calculational purposes. Instead, a more "friendly" geometrical version results if we 
perform a unitary transformation on ip defined by^ 

ip (t, x) = Uip (r, x) 
ip (r, x) = Ip (r, x) 

with the rotation 

u{R,{<p)Ry{e)) = u{R.,mu{Ry{e)) 



(2.3) 



exp 




exp 




(2.4) 



Geometrically, this amounts to a realignment of the coordinate axes with the chosen curvilin- 
ear coordinates (r, 6', 0). Under this field redefinition, the Euclideanized Lagrangian becomes 



Ip (r, x) Qip> (r, x) = ip (r, x) U'^QUip (r, x) = ip (r, x) Qip (r, x) 



(2.5) 



This transformation leaves the path integral measure invariant. 
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where 

e = 7^5o + 7l (dr + -]+ ih- (de + ^ cot ^^'j + ^l^^d^ + m , (2.6) 

and we used the fact that spatial rotations [in our case, U in Eq. (12. 4p ] commute with 7^;. 
Effectively, this choice leads to the selection of self-adjoint operators associated with the 
given coordinates. With our representation of the Dirac matrices (Appendix , we have 



e 





(2.7) 



sin 9 



-o"i / \ — 0"2 

It is now necessary to attempt a separation of variables to reduce the problem to one 
dimension. In particular, owing to the symmetry of the problem in angular variables, suggest 
an expansion of the form^ 

\Dn{r)Y2{9,4>) J 

Following Appendix |Bl we introduce the "generalized spherical harmonic" eigenfunctions Yi 
and Y2, 



if) oc 



(2.8) 



( 





1 




1 






de + - cot I 



-i 

1 





where 



A± = ± (j + 1/2) . 
Then, the effect of the Euclidean Dirac operator on the fields (12.81) is 



(2.9) 



(2.10) 



[ -^\{^r + l)C + ^D'^Y, + i-^tu^ + m)AY^ \ 

-{dr + }.)D + ^C] Y2 + [-lUJn + m] BY2 



Yi + (icOn + m) CYi 



(2.11) 



V ^[-{^r + \)B + ^A)Y2 + {^uJn + m)BY2 ) 



For notational simplicity we are omitting the j and mj labels on the spherical functions i^(i,2) in Eqs. ([2 
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We use the above harmonics to expand the spatial part of the fermion field, leaving the 
Euchdean time as before, 



3,mj,n 



j,mj,n 



-^nj(^)^2*fmj 



(2.12) 



Using 

II ^Ijmyij'm'. sin 0ded(f) = SjfSm^rn'. 

the Euclideanized action becomes 



Sf: — 



n,j,mj 



+ 



+ Cl^[dr + -\A^,-DlAdr + -]B^, 



'-nj ^nj 



+ {AnjAnj + B*jBnj — C*jCnj — D^jDnj) | . 



This expression can be written as 



(2.13) 



(2.14) 



n,j,mj 



(2.15) 



where 



njruj 



( \ 



B. 



Cnj 
V Dnj J 



(2.16) 



a 



njrrij 



nj nj nj nj 



—I 



+ im 
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leading, up to an irrelevant constant, to 



Z = JJ det 



(2i+i) 



-27^°) (w„) - 27(3) {dr + -] + Y^'— — — + m 



r 



We again use 



det (D) = Vdet (D) det (Dt) 
(up to a phase factor) and the properties of the Dirac matrices to get 



(2.17) 



(2.18) 



Z = JJ det 



(j + l/2) 



2^2 I d f ,d\ ^ {j + i/2y 
ar \ ar 



(j+1/2) 



^3 

as 



(2.19) 



Notice "det" in Eq. fl2.17p refers to both functional (r) and matrix sense. The operator in 
Eq. fl2.19p can be written as 

(j+1/2) 







72 ^3 

(j + 1/2) 

i^{rj) 0^3 



with 



Therefore, 



^ d f ,d\ (j + l/2f 



^{r,j) =^r, -— I r — I + 



det A 



" dr \ dr 

det^ 
det^ 



+ m 



(2.20) 



(2.21) 



Cz+i/2) ■ 











= det^ (A2 



(2.22) 



where 



+ ir,j) 



(j + 1/2) 



A, 



0'.'') 



(2.23) 



^+(r,j) 

As j takes only semi-integer values (j = 1/2, 3/2, 5/2, . . .), let / = j + 1/2 = 1, 2, 3 . . .. With 



this notation, 



1 d f 2d\ l{l±l) 



+ m . 



We can then write Eq. (12.191) as 

J] det ^'A^''"^) = Y[ {detQif J] (det 



s2l 



(2.24) 



(2.25) 



n,l 



n,l=l 



n,l=l 
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Now, 



H (det ni.,f = H (det nif^' = H (det ilif^' J] det fiz 

1=1 1=0 1=0 1=0 



Therefore, 



(2.26) 
(2.27) 



Z=(ll (detfi,)''^') ■ (2.28) 

\n,l=0 J 

As is well known, JJ^ (det QiY'''^^ is the spherical-coordinate version (up to a constant) 



n,l=0 



of detn^;, where De = d"^ + V^. In momentum space, Eq. fl2.28p then becomes, up to a 
constant, equal to the expression of Eq. (11. lip . Q.E.D. 



III. CONCLUSIONS 

We have shown the equivalence of the calculation for the partition function of a free gas 
of fermions in fiat spacetime in Cartesian and in spherical coordinates. The latter involved 
an expansion of the fermion fields in "generalized harmonics" as detailed above. While the 
result is not surprising, the treatment of the path integral in spherical coordinates in this 
paper presents novelties and subtleties hitherto not seen or emphasized. It is technically 
remarkable for instance, to see the emergence of the zero mode (/ = 0) for the scalar case in 
the process of "squaring" the fermionic determinant. Beyond its own merits, this calculation 
lends strong support to the correctness of the similar calculation in the case of a black hole 
background performed recently by the authors j4|. In that case, the result similar to Eq. 
f l2.19p in this paper after following the same procedure is 

2ff" - P 



J]det^+i 



/2 



+ + 2) / dr^ 



4 

0-3 
as 



dr \ dr I 



(Ti 

-ai 



(3.1) 



In Eq. (13. ip the function f = f[r) is defined by the Schwarzschild metric 
ds^ = f (r) ^ -g{r) dr^ - {dO^ + sin^ 6 d(j)^) , 

where 

r g[r) 



(3.2) 



(3.3) 



8 



in Minkowski spacetime with signature (+, — , — , — ). As we saw above, the determinant in 
Eq. fl2.19p has a very tight structure that allowed for a full diagonalization into products of 
determinants over subspaces of the original one. The presence of the last term in Eq. (13.1 p 
makes this process considerably more difficult than the flat spacetime case. The authors' 
NH expansion method made it possible to systematically calculate and isolate the leading 
contribution to the divergent part of the free energy — hence of the entropy — of the fermionic 
thermal atmosphere surrounding the black hole, which led to the Bekenstein-Hawking law 
within the framework of 't Hooft's brick- wall approach. Therefore, while the black- hole 
calculation is consistent with the expected fundamental result of black hole thermodynamics, 
it is important to establish the full validity of the technical aspects of the NH expansion, 
and this paper lends further credibility to our approach. 



Appendix A: Euclidean Dirac Matrices 

We define the Euclidean Dirac matrices to be 

= (Al) 

7| = -^7^ 

where, in our convention, the 7'^'s are the Dirac matrices for Minkowski spacetime with 
(+, — , — , — ). Of course, the signature of Euclidean space is (+, +, +, +), with the matrices 
satisfying 

{7^,7^} = 2<5^^ (A2) 
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Explicitly, in our representation, 



1% 



(l 





\ 




/o 








-^\ 


1 
















—i 








-1 










i 








\0 





-V 












/ 


( ° 





-l\ 




/o 





—i 


o\ 





1 
















i 





1 





7| = 


i 
















J 




Vo 


—i 





0/ 



(A3) 



Appendix B: Reduction in Spherical Coordinates 



Due to the structure of the Dirac operator, the separation into radial and angular com- 
ponents in this case is more complicated than in the one-component scalar case. One way 
to proceed is as follows: Eq. f l2.7p shows that the problem reduces to a two-dimensional one, 
so we postulate the "generalized" eigenvalue problem^ 



( 





1 




1 






+ i cot I 



-i 

1 





Solving for Y2 one can then obtain an equation for Yi: 



{de + -cot0)- . 

2 smtJ 



[dg + lcote) + -!—d^ 

2 smtl 



Writing 



yields the following 



m cos 6 



df n 1 

With the further change of variables, 

X = cos 6 

/ = (1 - X)^2-(1 + X)^2-P(x) . 



(Bl) 



(B2) 



(B3) 



(B4) 



(B5) 



Notice that A+ cannot be equal to A. 
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one obtains the following equation for P{x): 



X 



+ 



1 - {2m + 2)x 



dP 

dx 



+ 



A+A_ — mim + 1"; 



P = 0, 



(B6) 



This equation describes the Jacobi polynomials, for which A+A_ = — (j + 1/2)^ has to be a 
negative integer: A+A_ = —1, —2, —3, . . . (and hence j a semi-integer: j = 1/2, 3/2, 5/2, . . .). 
We can then choose X± = ±(j ' + 1/2). Solving for Y2 would give a similar equation to 
Eq. ( ]B6p but with the coefficient of the first derivative replaced now by —1 — (2m + 2)x. 
Therefore, we get 



Y ) 









X) 2 

m + 




(B7) 



In addition, these eigenfunctions y's are chosen to be orthonormal. 



II^2jmy2j'mr Sm 



^jj'^nijm'^ 



(B8) 
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